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Complete population transfer in 4-level system via Pythagorean triple coupling 
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We describe a relation between the requirement of complete population transfer in a four-mode 
system and the generating function of Pythagorean triples from number theory. We show that 
complete population transfer will occur if ratios between coupling coefficients exactly match one 
of the Pythagorean triples (a, b, c) C Z, c 2 = a 2 + b 2 . For a four- level ladder system this relation 
takes a simple form (V12, V23, V34) ~ (c,b,a), where coefficients Vij describe the coupling between 
modes. We find that the structure of the evolution operator and the period of complete population 
transfer are determined by two distinct frequencies. A combination of these frequencies provides a 
generalization of the two-mode Rabi frequency for a four-mode system. 

PACS numbers: 02.10.De, 02.20.Qs, 32.80.Qk, 42. 50. Hz 



Revealing hidden mathematical structures behind 
physical phenomena is of great importance, especially 
when qualitative dynamical properties of a quantum sys- 
tem reduce to a basic relation from the number the- 
ory. Here we show how the Pythagorean triple, which 
is the set of three integer numbers (a, b, c), satisfying the 



Pythagorean equation a 



c , is found to play a 



significant role in the dynamics of a four-mode system. 

Describing the evolution of a general multi-mode sys- 
tem and finding conditions for complete population 
transfer from one mode to another is a subject of exten- 
sive research for a variety of classical and quantum sys- 
tems. Coherent manipulation of population of states in 
atomic and molecular quantum systems [lj, spin control 
in nuclear magnetic resonances [2| , quantum information 
processing and directional optical waveguide tech- 

nology [5| are only a few examples where complete pop- 
ulation transfer is desired. Here we study the four-mode 
dynamics, which is of particular importance for the quan- 
tum information processing technology, where two-qubit 
quantum logic gates serve as elementary building blocks 
for designing fully functional scalable devices [H, |g, 0] • 

In general, solutions of dynamic coupled equations are 
difficult to analyze, and even for the simplest case of a 
two- level system, realized by a spin-i particle or a two- 
level atomic system, only a handful of analytical solutions 
are known in the literature. The simplest one, known as 
the Rabi solution describes a two-level system with 
a constant coupling. For this solution, complete popu- 
lation transfer between two states occurs only when the 
frequency of an external driving field is on resonance with 
the energy difference between the modes. The period 
of complete population transfer, called the Rabi flop- 
ping time, is inversely proportional to the strength of the 
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FIG. 1: Finite level realizations of the two- and four- mode 
systems of coupled equations (a) A four-level system with 
periodic nearest neighbor coupling, known as diamond shape 
structure, (b) A four-level system with V14 = 0, known as 
a ladder type structure, (c) a simple two-level system with 
coupling coefficient V eg and the detuning of A. 



mode coupling. A geometric visualization of a two-level 
system by representing its state by, a point within the 
Bloch sphere, plays important role in developing a clear 
intuitive understanding of the two-mode dynamics 0, [^] . 
Our approach helps to extend this geometric approach to 
four-level systems. 

The task of finding schemes for complete population 
transfer between selected states becomes increasingly dif- 
ficult in multi-mode coupled systems. There are group- 
theoretical methods which provide a rigorous tool of how 
to determine whether a system is wavefunction control- 
lable, i.e. when any initial state in a quantum system 
can be transferred into an arbitrary final state [loL [ill ] , 
however these methods are nonconstructive and do not 
provide a general recipe on how to implement complete 
population transfer scheme for a concrete system. So 
far, there is only a limited number of systematic meth- 
ods which can provide this goal. The schemes exploit- 
ing adiabatic evolution are known to be able to achieve 
this goal asymptotically, however these methods require 
strong pair-vise sequence of coupling pulses as well as 
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very long dynamical time [12, H3, H, [HI • Few solutions 
of the complete population transfer problem, which re- 
quires a set of coupling coefficients to satisfy some special 
relations, were found for N-mode systems [l6l [T7l] . 

In the present paper we describe a new analytical 
solution for a four coupled mode system with nearest- 
neighbor coupling. We show a clear similarity in the 
structure of this four-level solution and the structure 
of the Rabi two-level solution, originating from a com- 
mon geometric character of both solutions. We exploit 
the fact that a four-dimensional Hilbert space can be 
represented as a tensor product of two two-dimensional 
Hilbert spaces. For the four-level nearest-neighbor cou- 
pling system the relevance of this tensor-product struc- 
ture is revealed in its full simplicity when the Hamil- 
tonian is rewritten in the basis of the Bell states. To 
do a qualitative physical analysis of the resulting set of 
two-level equations we use a transformation known as 
the Hopf map [H, . The final equations for the evo- 
lution of two separate pairs of modes «-► \ipa) and 
\ip2) ^ | "04) provide natural generalization of the Rabi 
two-mode solution for thefour-level system. Finally we 
derive an equation for the time required for complete 
population transfer to occur. 

The requirement of complete population transfer in 
a four-level system imposes certain analytical relations 
on the coupling coefficients and we find that these rela- 
tions have special algebraic character: they are identical 
to a formula, which generates all primitive Pythagorean 
triples (PPTs). Such a triple is a set of three real num- 
bers a, b, c, which do not possess a common factor and 
satisfy the equation a 2 + b 2 = c 2 . For instance, (3; 
4; 5) and (5; 12; 13) are primitive triples, whereas (6; 
8; 10) is not PPT. Note also, that in spite of the fact 
that a set of numbers (1,1, V2) satisfy the Pythagorean 
relation, these numbers are not a Pythagorean triple. 
Through the centuries finding a formula to generate 
these triplets has intrigued both amateur and profes- 
sional mathematicians. The first general solution was 
given by Euclid in his Elements [2(|. It states that 
for any pair (p, q) of positive odd integers with p > q, 

is Pythagorean. 



the triple (a, b, c) = 



2 »P«> 2 

Other types of generating functions of PPTs can be found 
elsewhere [II Si, S3]- 

We start with a nearest-neighbor coupling four-mode 
Hamiltonian, when each level is coupled to two neighbor- 
ing levels, as in the diamond-shaped structure, shown in 

Fig. m 
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As an example of physical realization of this system, we 
can consider a laser-field driven four-level atom. Then in 
the rotating wave approximation the coupling coefficients 
are defined as V%j — i-iijC (t) /h. Here e(t) is the field 



amplitude and /i^- are dipole matrix elements between 
nearest- neighbor levels i and j. A general case of periodic 
nearest- neighbor coupling (diamond-type structure), is 
schematically shown in Fig. [2a), and particular case of 
V\4 = 0, describing a ladder-type coupling, is presented 
in Fig. 

We rewrite Hamiltonian (fTJ) in the Bell basis using the 
unitary transformation 
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The wave-functions are transformed as \ip B 



nJ - W\1> n ), 

and the Hamiltonian becomes Hw — W^HW. In this 
basis two distinct su(2) subalgebras, formally associated 
with qubits 1 and 2, can be identified. We rewrite the 
Hamiltonian as a linear combination of Pauli matrices 
&(x,y,z) acting on qubits 1 and 2 



H w = ftd) ® /(2) + /(i) ® ft(2) , 

1(1) _ (V23+V14) - 1 (Vl2-V 3 4) ~ 

11 — 2 u z \ 2 x ' 

1,(2) _ (V12+V34) - (V23-V14) ±. 



(3) 



In the language of Lie group theory Hamiltonian (|TJ) 
generates a subgroup SU{2)®SU{2)/Z 2 of the full SU{4) 
group. The dynamic problem factorizes into two sep- 
arate problems for two SU(2) unitary operators acting 
on two qubits, thereby geometric tools for visualization 
of resulting solutions are readily available. By using an 
algebraic property of local transformations, we can rep- 
resent the action of an SU(2) x SU(2)/Z2 operator on a 
four-dimensional state vector as left and right multipli- 
cation by two 2x2 SU(2) matrices, acting on a 2 x 2 
complex matrix. The latter represents an element of the 
four-dimensional Hilbert space [HI • Thus we rewrite the 
equations for the evolution of amplitudes a n (t) of the 
states \ip n ), n £ {1, 2, 3, 4} in a form of two rotations: u% 
acts from the left, and u 2 acts from the right: 



A (t) = ux (t) A(t = Q) ul (t) 



(4) 



Here A (t) contains the information about four ampli- 
tudes 01,2,3,4, 

A (t) = ax (t) I + a 2 (t) a x + ia 3 (t) a y + a 4 (t) a z , (5) 

and the operators ux and u 2 are local rotations of qubits 
1 and 2 correspondingly 



;/ , (/) = exp { i l - [(Via - + ( v ™ + °* 



u 2 {t) = exp\i-[{V l2 + V ?>A )d x -{V 23 -VxA)o- z \ 



It immediately follows from this relation that if the 
system is initialized in the ground state |'0i) so that 
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TABLE I: Comparison of the Rabi solution for two mode system (the middle column) with a four-mode nearest neighbor 
solution (the right column) . We observe A striking similarity in the structure of both solutions. 



Parameter 


Two-mode dynamics 


Nearest-neighbor four mode dynamics 


Dynamical symmetry 
Spanned space 

Generalized frequencies 

"Torque" vector 
Ground state evolution 

Excited state evolution 

Inversion time 


SU{2) 

i\i> B )Me)} 

V = y/V& + 

tl R = (Re{V 12 },Im{V 12 },A) 
a g (t) = cos {Vt) - y==^f sin {Vt) 

T = T^TT 


SU(2) x SU{2) 
{|^>,|V3>} and{|^ 2 ),|V4>} 

V L = y/(Vv, ~ V 34 ) 2 + (V 23 + V 14 ) 2 
Vr = sj (V u + V34) 2 + (V23 - V14) 2 

fip = ^ (6,6,6) 

01 it) - cos (V L t) cos (V R t) - | 3 2 sin (V L t) sin (V R t) 
a 3 (t) J£= sin (V L t) sin {Vat) 

V s i +'3 



ai (0) = 1 and (12,3,4 (0) = 0, the amplitudes of ai ; 3 (t) 
will remain real and the amplitudes of 02,4 (t) will be 
purely imaginary. 

To calculate explicitly a n (t) and analyze the problem 
of complete population transfer between the states 
and I "03), we first note that there is another symmetry, 
which allows us to solve the problem in an elegant geo- 
metric fashion. Suppose that we chose to rotate the basis 
in the space spanned by vectors \ip2) and \ipi). Such rota- 
tion apparently will have no effect on the evolution of the 
states l^i) and |^is). We can represent this transforma- 
tion as a phase multiplication acting on two complex vec- 
tors (Via + iVu) (Via + iV l4 ) e lf> and (V 23 + iV 3A ) 
(V23 + iV^) e %B . The invariance of the amplitudes a\ (t) 
and 03 (i) under such a transformation means that the 
amplitudes a\ and 03 are determined not by the full set 
of coupling coefficients {V12, V23, V34, V14} € i? 4 , but by 
an element of the quotient space R 2 x R 2 / SO{2), de- 
scribed by a special algebraic transformation, known as 
the Hopf S 3 — ► S 2 projective map [IHIIH . In physics, the 
Hopf map is commonly associated with the Bloch Sphere 
representation of a pure state. In our problem, the map 
takes the 4-dimensional Vy space to a 3-dimensional 
space, 
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V12V23 + V14V34 , 
V12V34 — V23V14 , 

\ (v 2 2 + V 2 , - v 2 3 ~ v 3 \) 



(0) 



The new coordinates satisfy the equation £q — £? — 
£2 — £3 = 0, which is the equation for a 3-dimensional 
cone embedded in four-dimensional Eucledian space [25| . 
Physical meaning of this set of parameters will be clari- 
fied when we use them to derive the final solution for the 
amplitudes a n (t). 

Time dependent amplitudes of the modes \ipi) and ^3) 
can be expressed directly from equation ([4]) using cou- 
pling coefficients Vjj , however by using the projective 
coordinates {^1,^2,^3}, algebraic expressions for these 
amplitudes can be significantly compactified. As it im- 
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FIG. 2: Complete population transfer ji/>i) < > \tjj3) (or 

I ^2) < > l^})- (I) Level populations as functions of inter- 
action time, for a ladder-type coupling scheme. The system 
initially is prepared in the ground state (solid red) . The pop- 
ulation is periodically transferred to the third state (dashed 
green). Parameters correspond to {p, q) = (3, 1) in Eq. 8. (II) 
The Pythagorean triple relation between the coupling coeffi- 
cients in four level ladder system. The computed transition 
time for complete population transfer matches the calculated 
transition time from Eq. 9. 



mediately follows from equation (|4|), 

: sin {V L t) sin (V fl t) , 



£1 



cos (V L t) cos (V R i) - f sin [V L t) sin (V R t) . 
V Si + ?3 



a 3 (i) 
01 (f) 



Here Vl,r = y 5 (£0 T £2) are the left and right frequen- 
cies, respectively. Explicit formulas for Vt and Vr are 
given in Table I. 

The form of the exact solution given by equation (|7|) 
is similar to the form or the Rabi solution for a two-level 
system, as demonstrated in Table I. In the two-level case, 
shown in Fig. 1(c), we use notations a g and a e for the am- 
plitudes of the ground and excited states, while detuning 
was denoted as A = {ui e — ui g ) — uji aser . As can be seen 
from Table I, instead of one frequency V for the two- 
mode dynamics, the four-mode system is characterized 
by two generalized frequencies Vl and Vr, and instead 
of single sine and cosine functions in the two mode case, 



(7) 
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we see that the four-mode dynamics is determined by the 
product of two sine and two cosine functions. The pre- 
factors of the sine function in both cases share a similar 
structure, being formally equal to x- and z- components 
of unit vectors. One can see that there is a very close 
analogy between variables £1,2,3 and components of the 
the torque vector for two-mode dynamics 

Now we have all the necessary equations to solve the 
problem of complete population transfer. To realize com- 
plete population transfer from the state to the state 
^3) at time t = t one has to set the £3 variable equal 
to zero, in the same fashion as the detuning A in the 
two-mode case has to equal zero (the requirement of on 
resonant interaction). Next, complete population trans- 
fer will occur only when dynamic angles Vti and Vpt, 
simultaneously complete a 7r-phase rotation, i.e. when 
V L = § ? (2m! + 1) = £p and V R = £ (2m 2 + 1) = £q. 
After some trivial algebra, we derive the following solu- 
tion: 

ir 2 I p 2 + q 2 p 2 + q 2 \ u 2 

(£o,£i,6) = ( — 2 — ,pq ' — 2 — J = T ( c ' a ' fc ) -( 8 ) 

This solution exactly matches the definition of t he g ener- 
ating function of primitive Pythagorean triples [201 ] . For 
the nearest-neighbor 4-mode coupling problem complete 
population transfer can occurs between two nonadjacent 
states I -0i } and |V>3) (or \1jj2} and IV^)) only when the ratio 
(£0 : £1 : £2) is equal to ratio of a Pythagorean triple. The 
equation for coefficients {V12, V23, V34, V14}, can be ob- 
tained from equation (8) by inverting the transformation 
(6). For the special case of latter-type coupling, where 
V14 — 0, this solution becomes £0 = V 2 2 , £1 = V23V12, 
£2 = V34V12 such that relation between the nearest- 
neighbor coupling coefficients and the Pythagorean triple 
takes a simple form of a proportion (V12; V23; V34) ~ 
(c,a,b). 

We tested our theoretical prediction by performing nu- 
merical simulations on the dynamics of a four-level ladder 



transitions in i?&85 : 5<Si/2 «-» 5P3/2 «-> 4_D 3 / 2 4^/2, 
with resonant CW interaction of 780.2r7.rn, 1.529/jm, 
1.344/im, respectively. The coupling coefficients were 
chosen to satisfy the simplest Pythagorean triple ratio 
(V12 : V23 : V34) ~ (5 : 3 : 4). As seen in Fig. 2, numeri- 
cal results are in complete agreement with the analytical 
solution and confirm that there is periodic population 
transfer between states and ^3). The time period 
for complete population transfer is given by 



sip ^/vl+vi V*To ■ 

Here we denoted the transition time required to achieve 
population transfer as Sip, analogous with the Rabi fre- 
quency for a two-level dynamics. Note that this param- 
eter scales as the absolute value of the torque vector, 
similar with the two-mode case. 

In conclusion, we identified a new scheme for complete 
population transfer in a four-mode systems. We observed 
very close connection between the structure of solution 
for the nearest-neighbor coupling four-level system and 
the generating function of primitive Pythagorean triples. 
This solution can be used not only for time-dependent 
problems, but also for some problems of spatial propaga- 
tion of light pulses, such as coupling between directional 
waveguides. We expect that similar solutions, revealing 
deeper link with the number theory, can be found for 6- 
and 8-level systems. The present method, describing the 
four-level dynamics, can be generalized to include more 
complex exactly solvable two-level models. This work is 
in progress. 
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